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Abstract 

A new method is proposed for practical calculation of the effective interaction between impurity scatterers in superconductors, 
based on algebraic properties of related Nambu matrices for Green functions. In particular, we show that the density of states 
within the s-wave gap can have a non-zero contribution (impossible either in Born and in T-matrix approximation) from 
non-magnetic impurities with concentration c< 1, beginning from ~ c 3 order. 



I. INTRODUCTION 

Impurity effects are at the center of interest in stud- 
ies of superconducting (SC) materials, especially of those 
with high transition temperature (HTSC). In the gen- 
eral theory of disordered systems with disorder due to 
diluted impurity centers, the so-called group expansion 
(GE) method was proposed as most consistent for quasi- 
particle Green function (GF) Q, and it was also formu- 
lated for SC systems Q. 

The GE's of different types (see below) are analogous 
to the classical Ursell-Mayer group series in the theory 
of non-ideal gases 0], where the particular terms (the 
group integrals) include physical interactions in groups 
of the given number of particles. In the quantum theory 
of solids, GE includes indirect interactions (dependent 
on the excitation energy e) between the impurity cen- 
ters, through the exchange by virtual excitations from 
(admittedly renormalized) band spectrum, so that each 
term corresponds to summation of certain infinite series 
of diagrams. These expansions were elaborated in a detail 
for various kinds of normal quasiparticle spectra, where 
they define the interplay between extended and localized 
states p|, however their usage in SC systems encoun- 
ters considerable technical difficulties due to existence of 
anomalous GF's. 

The present paper is aimed on an efficient algorithm for 
resolving these difficulties. We develop the specific alge- 
braic techniques to calculate Nambu matrices in various 
terms of GE for the exemplary case of s-wave symmetry 
of SC order parameter, permitting to explore the impu- 
rity effects in this case beyond the scope of Anderson 
theorem In particular, we find that pair clusters of 
impurities (2nd term of GE) can not produce finite con- 
tribution into the quasiparticle density of states (DOS) 
within the s-wave gap, alike the simplest effect of iso- 
lated impurities (1st GE term), but a non-zero contribu- 
tion into the in-gap DOS is already possible for the 3rd 
GE term (impurity triples). These results allow also a 
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straightforward extension to the d-wave symmetry char- 
acteristic for doped HTSC materials (where the dopants 
not only supply the charge carriers but also play the role 
of impurity scatterers) . 



II. HAMILTONIAN AND GREEN FUNCTIONS 

For description of electronic spectra in a SC system 
with impurities, it is convenient to use the formalism of 
Nambu spinors: the row-spinor ^ = (a k *, a-k,j.) and 
respective column-spinor ifa, writing the Hamiltonian in 
a spinor form 

k p.k' 

(1) 

It includes the normal quasiparticle energy £k, the mean- 
field gap function Ak, the Pauli matrices Tj and the 
perturbation matrix V — Vi/tq. The impurity (attrac- 
tive) perturbation on random sites p with concentration 
c = ^ p 1/N 1 is described by the Lifshitz parameter 
V L . 

The energy spectrum of a Fermi system is described 
through the Fourier transformed two-time Green func- 
tions (GF's) 0: 

({a\b)) e = i f e^- iO »({a(t),b(0)})dt, (2) 

J — oo 

where (. . .) is the quantum statistical average and {., .} 
is the anticommutator of operators in Heisenberg repre- 
sentation. For the system, Eq. ^ we define the 2x2 
Nambu matrix of GF's 

G k ,k' = ((V^k|<». (3) 

The matrix elements in the expanded form of Eq. are 
the well-known Gor'kov normal and anomalous functions 
Q. In what follows, we shall also distinguish between 
the Nambu indices (N-indices) and the quasi-momentum 
indices (m-indices) in this matrix. In absence of im- 
purities, the explicit form of the matrix, Eq. 03 turns: 



Gk,k' — * <5k,k'G k , where the non-perturbed (m-diagonal) 
GF matrix 



G k = 



£ + 6c t 3 + A k ?i 



El 



(4) 



involves the SC quasiparticle energy i? k = + A k . 
The relevant physical properties of SC state are suitably 
expressed in terms of these GF's. For instance, the single- 
particle DOS, related to the electronic specific heat, is 
given by 



k 



(5) 



where Gk = G kjk is the m-diagonal GF. 

Now we pass to calculation of GF's in SC systems at 
finite concentration c of impurity centers and analyze ex- 
plicit structure of corresponding GE's. 



III. GROUP EXPANSIONS FOR SELF-ENERGY 

We derive GE's for the system defined by the Hamilto- 
nian Eq. starting from the Dyson equation of motion 
for a matrix GF: 



Gk,k' = G k £k,k' 



e i(k "- k> ^re k «,k'. (6) 



P ,k" 



A routine consists in consecutive iterations of the same 
equations for the GF's in the "scattering" terms of Eq. 
Eland separating systematically those already present in 
the previous iterations jlf]. Thus, for the m-diagonal GF 
Gk, we first separate the scattering term with the func- 
tion G k itself from those with G k ' k , k' ^ k: 

G k = G k - 1 e I ( k - k ') p G k i>G k ', k = 



k', P 



= G k -cG k i>G k -l J2 ci(k_k ' ) p Gk^G k ',k(7) 
k'#k, P 

Then for each Gk',k, k' ^ k we write down Eq. again 
and single out the scattering terms with Gk and Gk',k in 
its r.h.s: 



Gk'.k — — 



-cG 



.1 Y, e l ( k '- k ") PG k ,yG k ^ : 



k", P ' 



N 



-1 y e < k '- k ) p 'G k ,yG k 



p'#p 



.1 e< k '- k ")-PG k ^G k » k . (8) 



k"^k,k';p' 



Note that among the terms with Gk, the p' = p term 
(the second in r.h.s. of Eq. |HJ bears the phase factor 



e i(k k) P ^ g0 ^. j g co jj eren t to that already figured in the 
last sum in Eq. [3 That is why this term is explicitly sep- 
arated from other, incoherent ones, oc 
(but there will be no such separation when doing 1st it- 
eration of Eq. Elfor the m-non-diagonal GF Gk",k itself). 

Continuing the sequence, we collect the terms with the 
initial function Gk which result from: i) all multiple scat- 
terings on the same site p and ii) such processes on the 
same pair of sites p and p' ^ p. Then summation in p of 
the i)-terms gives rise to the first term of GE, and, if the 
pair processes were neglected, it would coincide with the 
well known result of self-consistent T-matrix approxima- 
tion @. The second term of GE, obtained by summation 
in p,p' 7^ p of the ii)-terms, contains certain interaction 
matrices ^4 P ', P generated by the multiply scattered func- 
tions Gk'.k, k' ^ k, etc., (including their own renormal- 
ization). For instance, the iterated equation of motion 
for a function Gk",k with k" ^ k, k' in the last term of 
Eq. will produce: 



6k»,k = ~ e< k "- k "')' p "G k ,,\>G k ^ 



N 

k"', P " 

e »(k"-k)-p e i(k"-k). P ' 



k — 



-G° k „VG k 



Gl„VG k 



N N 
+ terms with G k ' k and G k " k + 
+ terms with G k », k (f/kk' k"). (9) 

Consequently, we obtain the solution for an m-diagonal 
GF as 



Gk — Gk,k — 



(at)" 1 



(10) 



with the matrix GE for the renormalized self-energy ma- 
trix: 



£ k = cT 



l + c> U), ln e 

X (l - A),n-4n,0 



ikn 



(11) 



Here T = — V^l + GVj is the renormalized T- 

matrix, and the indirect interaction (mediated by the 
quasiparticles of host crystal) between two impuri- 
ties at lattice sites and n is described by the ma- 
trix l ,n = N- 1 Ek'/k eik n Gk'f , with the sum in 
quasi-momenta restricted due to the above algorithm 
of separation. There are even more such restric- 
tions in each product of these matrices: Ao jn ^n,o = 
^- 2 E k ^kEk^k,k'e^ k '- k ") n Gk'TG k ,f, and so on. 
This seems to seriously hamper calculation of the sum 
En^o m Eq- El (not to say about higher GE terms). 
However, the difficulty is avoided, taking into account 



2 



the identities for two first terms of its expansion Q : 



^0 + ^EE el(k '^ k)n ^T: 



n k'^k 



-A 



0,0; 



and 



^0,n^4n,0 — — j4q,0 + E A),n-An,0 — —^0,0 



n k',k"#k' 



- -A 2 



due to the momentum independence of T-matrix, and 
the fact that the restrictions can be simply ignored in 
the higher products, like 



ikn 



-A 3 



n^0 



+ Y ^0,nA n ,oA),ne lkn = -A 3 Q + 
n 

^E E e ! ( k '- k " +k "'- k ) n Gk'fGk»fG k -f = 



n k',k"^k' 



r2 E ^k'^Ck'-rGk-k'+k"? 1 , 



^0,0 ■ ^2 

k',k" 

etc. Thus we arrive at the final form for the renormalized 
GE 



S k = cT [l - &A ,o - cAl + 

+^E(^n^' 



1 — ^0,n 



12) 



where Aq ^ = Gq :II T and the renormalized local GF ma- 
trices G , n = N^ 1 X] k eJkn Gk and G = G ,o are already 
free from restrictions. The two terms, next to unity in 
the brackets in Eq. 1121 correspond to the excluded dou- 
ble occupancy of the same site by impurities, the sum in 
n^O describes the averaged contribution of all possible 
impurity pairs, and the dropped terms are for triples and 
more of impurities. 

An alternative routine for Eq. ^consists in its iteration 
for all the terms Gk",k and summing the contributions 
oc G k , like the first term in r.h.s. of Eq. El This finally 
leads to the solution of form 



fiO i /~i0 v*o fiO 

k — <-fk ' "k^k^ki 

with the non-renormalized self-energy matrix 



(13) 



E k = cf 




All 



!-(<„ 



(14) 



and the respective elements T° = — V ( 1 + G°V 



Ao, n — Gq^T , Gq — N 1 J2ic e lkn G k and G° — Gq Q . 

Presenting GF's in the disordered system in the form of 
GE's generally leads to respective expansions for its ob- 
servable characteristics. For instance, the impurity per- 
turbed DOS is expected in the form: p (e) = po (e) + 
pi (e) + p2 (e) + • ■ ■, related to contributions of pure crys- 
tal, isolated impurities, impurity pairs, etc. 

However, usage of each type of GE, the renormalized 
Eq. El or the non-renormalized Eq. ^] is only justified 
if they are convergent (at least, asymptotically). Since 
the matrices T and A are energy dependent, convergence 
of each type of GE is restricted to certain energy ranges, 
and these ranges are generally different. For a number 
of normal systems with impurities, where GE's are con- 
structed of scalar functions A) jn , it was shown that the 
renormalized GE converges within the region of band- 
like states, well characterized by the wave- vector, and the 
non-renormalized GE does within the region of localized 
states Q. To get quantitative estimates of convergence 
and higher order contributions to self-energy, operating 
with the matrix functions ^4o,n in Eqs. Eland^J a spe- 
cial techniques is necessary that we construct below. 



IV. ALGEBRAIC TECHNIQUES FOR NAMBU 
MATRICES 

Let us explicitly calculate the elements of above defined 
GE's for the simplest s-wave symmetry of the SC gap 
function: A s (k) = A. The unperturbed local GF matrix 
is obtained as an expansion in Pauli matrices: 



G y 



(e + An) g - g 3 T3, 



(15) 



where 



9o (e) 



irp-p 



defines the s-wave DOS in pure crystal: po (s) = 
(2/7r)Imgo (e) (with the Fermi DOS p-p of normal quasi- 
particles), and the electron-hole asymmetry factor 



53 [£) 



1 



is almost constant and real. 

Then we readily calculate the non-renormalized T- 
matrix: 



T° = 



T 3 + V- 



e- A?i 



np F l + v' 2 \'" ' " \/A 2 - e 2 , 
with the dimensionless perturbation parameter 
7T VlPf 



(16) 



v 



2 1 - Visa 



3 



Since the denominator 1 + v 2 in Eq. El can n °t be zero, 
the quasiparticle localization on a single impurity center 
in the considered s-wave superconductor turns impossible 
0. If the self-energy is approximated by the first term 
of GE: E k w cT°, then Eq. Housed in Eq. [Hand in Eq. 
pleads to the same DOS po (e) (that is, pi (e) = 0) with 
the same gap value A as in pure crystal. This justifies 
Anderson's theorem p| within T-matrix approximation 
for an s-wave SC with point-like impurity perturbations. 

However, even if there is no in-gap poles in the single- 
impurity T-matrix term, they can appear in the following 
terms of GE, that would describe localized states on im- 
purity clusters 0. Thus, the pair contribution: 



and then Eq. |22 implies an important formula for their 
arbitrary product: 



(24) 




M COS ^ Vn, , Sin ^ ^n; 

\ i=l i=l / 

This formula permits to reduce an arbitrary polynomial 
of A),n matrices to a single M-matrix whose arguments 
are polynomials of J-"o,n functions. The next important 
property of M-matrices is that the determinant: 



Pi (e) 



nN 



J2 Tr ImG°T n cos kn + (l° n ) ' 



k,n^0 



det 



1 - M (a, b) = (1 - a) 2 + b 2 



(25) 



can be zero only if the components are a — 1 and 6 = 
simultaneously. 



(17) 



should only follow from the poles of [1 — (^o in ) 2 ] 1 ; since 
the matrix Ag n (e) is real at e 2 < A 2 . 

In fact, the long distance asymptotics of this matrix 
(at n~^> a) is: 

^o,n ~ -^b,n (e) (cos ip n + £ ^ 3 A = ?A ^ 2 sin ip n J , (18) 



VA 2 - e 2 



where cp n = kp\n\ + 5 and the particular forms for the 
scalar "envelop" function !Fo, n (e) and the phase shift S 
depend on the system dimensionality: 



Fa,n (e) = 
for 2D, and 
Fo,n (e) 



—, cot<5=- , 19 

\/rr^ 2 v/tt/cf | n | i+« 



\ n \/ r ? 



VTTv 2 k ¥ |n| 



-, cot<J = l/u, (20) 



for 3D, wit h the en ergy dependent decay length r 6 = 
a 2 fc F /(7r jO0 VA 2 - e 2 ). 

The following analysis is essentially simplified, intro- 
ducing matrices of the structure: 



M (a,b) = a + b 



£?3 — tA?2 

VA 2 - e 2 ' 



(21) 



which form a closed algebra with the product rule for the 
a, b components: 

M(a,b)M{a',b') = M (aa! - W , ab' + bo!) . (22) 

In this notation, the interaction matrices, Eq. El are 
presented as 



^0 n ~ Fq,tiM (COS (p n , SHI (f n ) , 



(23) 



V. IMPURITY CLUSTERS AND IN-GAP DEN- 
SITY OF STATES 

The above developed techniques permit to quantify 
the effects of impurity clusters in quasiparticle spectrum. 
Thus, we conclude that the necessary condition for the 
pair contribution, Eq. ^3 



det 



det 



0(26) 



and J-Q n = 1. 



is only possible if (p n = irq, q = 1, 2, . 
But, taking into account the exponential factor e~' n '/ re < 
1 in Eqs. 1191201 this requires that: 

o „,2 



ttI 



> Trq — arctan 



1 + v 



1 



for 2D, or 



> Trg — arctan v, 



for 3D, which can not be fulfilled at any v and q > 1. 
Hence there is no contribution to DOS within the s-wave 
gap from impurity pairs (p2 (e) = 0), the same as from 
the single-impurity T-matrix, Eq. El 

But for the next, triple term of GEjwhich contains the 
inverse of the matrix (see, e.g., Ref. Q): 



1 — 2v4g n ylQ n ,A° n , — ^o : „) (^0,n') (^n,n') — 

= 1-M(a n , n ,,6 n , n 0, (27) 



the conditions by Eq. |2£lturn already realizable. They 
are explicitly formulated as follows: 

Z^O.n-^O.n'-^n.n' COS (i/3 n + (f n > + (£>„_„') + 
+^"o,n cos 2<fi n + Tl^ COS 2f n , + 
+^ 2 n ,cos2( /3n _ n , = 1, (28) 



4 



eral form for p^ (e) as an integral 
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Figure 1 : An example of simultaneous solution of the two con- 
ditions, a n n / = 1 and b n Y1 ' = 0, necessary for impurity triples 
to contribute into the in-gap DOS of a planar SC system. In- 
set: the chosen geometry of equilateral triangle 0, n, n'. 

and 

K.n' = 2:F , n :Fo,n':Fn,n' shl (p n + (f n > + ^„_„') + 

+^"o n sin 2tp n + Ti n, sin 2ip n > + 

+^ n , sin2 ^ n _ n , =0. (29) 

The easiest localization of quasiparticles of course is ex- 
pected at the very edge of the gap: e — *■ A, where 
e -l n l/»"e i Then, e.g., for v w 1.728 it is achieved 
with |n| = |n'| = |n — n'| sa O.SGGkp 1 , as shown in Fig. 
1, and it can be also reached for close values of v by small 
adjustments of the distances. 

Once such a pole exists, we can introduce in the con- 
figuration space of triangles 0, n, n', in the vicinity of this 
pole, the effective variables: 



r = J (1 - a a<n >) + b 2 nn ,, tp = arctan J 1 '" — , 
and a certain z, independent of r, tp, arriving at the gen- 



Pa (e) = c 3 Im / — L-li^ — ld r dp dz 
.1 r 

= 2ttc 3 J F (0, ip, z) dip dz. 

Here the function 

F (r, ip, z) = 2{a z , v cos <p + b z , v sin <p)J(r, p, z) 

contains the M-algebra coefficients a Z)lp , b Z)V related to 
a certain matrix numerator in the triple GE term (see 
its explicit form in Ref. and the Jacobian J(r, tp, z) 
of transition from n, n' to the effective variables. Apart 
from technical details of calculation, this contribution ~ 
c 3 pp can be routinely obtained. 

It should be noted however that the above approach 
uses the asymptotical form, Eq. 1181 of interaction func- 
tions, so it is only justified if the resulting poles, like that 
in Fig. m are related to long distances n 3> a. This is 
easier achieved if afcp -C 1 for the host crystal, as is the 
case for HTSC systems, but hardly for common s-wave 
superconductors. In the latter case, the finite in-gap DOS 
should result more probably from high enough terms of 
GE. 

Another important aspect is the analysis of conver- 
gency of matrix GE's, Eqs. 1121 1141 permitting to distin- 
guish between band-like and localized spectrum areas in 
the non-uniform SC system. This can be realized, using 
the above^estimates for the non-renormalized functions 
n and T° as approximations for the renormalized ones 

Ao, n and T in the region of band-like states e > A. 

A similar algebraic techniques can be also developed 
for the d-wave SC systems, related in that case to the 
impurity resonance states. In this situation, the higher 
DOS corrections p%,... to non-zero po + pi within the gap 
can be not so pronounced as for the zero-background case 
of s-wave system. Nevertheless the important criteria for 
GE convergence can be obtained, permitting a more con- 
sistent validation of T-matrix approximation and demar- 
cation between band-like and localized states, compared 
to the recently suggested check based on the phenomeno- 
logical Ioffe-Regel criterion 
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